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We study the non linear response of current-transport in a superconducting diffusive nanowire 
between normal reservoirs. We demonstrate theoretically and experimentally the existence of two 
different superconducting states appearing when the wire is driven out of equilibrium by an applied 
bias, called global and bimodal superconducting state. The different states are identified by using 
two probe measurements of the wire, and measurements of the local density of states with tunneling 
probes. The analysis is performed within the framework of the quasiclassical kinetic equations for 
diffusive superconductors. 



I. INTRODUCTION 



Superconducting nanowires are often part of objects 
to study the Josephson-effects in graphene, carbon nan- 
otubes or semiconducting nanowires. In addition in many 
cases superconducting nanowires themselves are used to 
study their response to radiation. In most cases the elec- 
tron back-scattering resistance is assumed to be located 
at the interfaces and in the normal metal part. An in- 
teresting question is to what extent the superconducting 
mesoscopic or (nano-)wires themselves contribute to the 
resistance of a device due to the conversion from nor- 
mal current to supercurrent and vice versa. For super- 
conducting nanowires between superconducting contacts 
a common assumption is that the applied power leads 
to dissipation and to an increased temperature varying 
over the wire length!^ In quite a few experiments with 
a nanowire between normal or superconducting pads, a 
parabolic temperature profile T(x) is assumed to con- 
trol the local superconducting propertiesPEl The defini- 
tion of a temperature however requires that the electrons 
are in local equilibrium, a condition not easily met for 
wires of mesoscopic length scales. In the case of a bi- 
ased normal wire^ the diffusion time, Trj — L 2 /D, with 
L the wire length and D the diffusion constant, is much 
shorter than the relaxation time for electron-electron in- 
teraction, r ee . In this case the electron distribution is 
highly non-thermal and given by a two step function 
f(E,x) = (1 - x)f {E - eV/2) + xf (E + eV/2), with 
fo(E,T) = l/(exp(E/kT) + 1) a Fermi-Dirac distribu- 
tion, eV the applied bias, k Boltzmann's constant, E the 
energy of the electrons measured from the Fermi energy, 
T the bath temperature and x the coordinate along the 
wire. A general, non-thermal (or non-equilibrium) elec- 
tron distribution in a superconductor influences almost 
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all aspects of that superconductor. It affects the local 
Cooper pair density and the current-carrying capacity, 
but it can also produce a voltage-drop in the supercon- 
ductor, i.e. a DC resistance of the superconductor. To 
discuss the various contributions it is advantageous to 
separate the non-equilibrium distribution function, f(E), 
into an energy (or longitudinal) mode, acting primar- 
ily on the amplitude of the superconducting gap, and a 
charge (or transverse) mode, /t, which lead to a shift 
in the pair chemical potential /z cp PThe latter mode fx 
describes an imbalance between electron and holes in the 
excitation-spectrum, leading to a net charge Q* in the 
(decaying) excitations. This contribution can be domi- 
nant in experiments probing electrical transport in su- 
perconducting heterostructures at subpgap energies. 
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FIG. 1: A superconducting Al nanowire connected to two 
massive normal reservoirs, consisting of the same Al, cov- 
ered by a normal metal Cu layer: (a) SEM-picture, (c) AFM- 
picture, (d, e) schematic representation. The thin Al of the 
pads is driven normal by the inverse proximity effect of the 
thick normal Cu. Normal tunneling probes are attached for 
local measurements (b). 
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In this paper we report on an experimental and theo- 
retical study of nonli nea r electrical transport in a well- 
defined model systemP^in which a superconducting wire 
is connected to two large normal contact pads (Fig. [T]). 
The normal electrodes induce evanescent subgap states in 
the superconducting wire. In addition they act as equi- 
librium electron reservoirs to fill and empty the states 
in the superconducting wire. When a bias eV is applied, 
evanescent electrons and holes are injected from the reser- 
voirs into the superconducting wire, and the resulting 
non-equilibrium distribution function consists of both an 
energy mode /l and a charge mode /t- The well defined 
boundary conditions and simplicity of this system make 
it a natural choice to study the superconducting state in 
the presence of a general non-equilibrium. 

We address these microscopic properties of the wire 
experimentally using two point measurements of the 
nanowire, which are a sensitive probe for the resis- 
tive properties originating in fx- Measurements with 
tunneling probes allow to measure the local density of 
states and the different chemical potentials involved. We 
demonstrate that two distinct metastable superconduct- 
ing states exist when the wire is driven. The first su- 
perconducting state extends over the complete length of 
the wire, and has been reported in the linear regime by 
Boogaard et alP The second state exists only under driv- 
ing, and consists of two, geometrically separated super- 
conducting domains, both at the ends of the wire. We 
show that the superconductivity nucleates in the vicinity 
of the normal reservoirs, because the local electron dis- 
tribution is closer to the equilibrium state. We analyze 
these states using the quasiclassical Green's functions, 
and show how the energy mode controls the existence of 
these states, whereas the charge mode controls the resis- 
tance. Hence, the full non-linear response is found to be 
the result of a complex interplay between both the charge 
and the energy mode non-equilibrium. 



II. THEORETICAL FRAMEWORK 

We consider a model system consisting of a super- 
conducting one dimensional diffusive wire connected to 
two normal, equilibrium reservoirs (Fig. [TJi). Electrons 
are injected into and extracted from the superconduct- 
ing wire by the reservoirs with equilibrium Fermi distri- 
butions fo(E ± eV, Tf,), with relative Fermi-levels deter- 
mined by the applied voltage V and Tf, the bath tem- 
perature of the reservoirs. Within the wire the elec- 
trons are distributed over the energies with a position 
and energy dependent non-equilibrium distribution func- 
tion f(E,x) determined by a diffusion equation. In addi- 
tion the electronic states concerned are decaying states, 
evanescent modes, as their energy is smaller than the en- 
ergy gap (eV < 2 A). Therefore it is necessary to include 
the interplay between these short-lived states and the 
superconducting condensate, which goes beyond a two- 
fluid description, in which a sharp distinction between 



long-living quasiparticle states and the condensate is as- 
sumed. Such an analysis is performed using the quasi- 
classical Green's functions theory for superconductivity, 
which treats the electronic properties of the excitations 
and the condensate on the same footing!^! 
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The retarded (advanced) functions G R ^ consist of nor- 
mal and anomalous propagators G and F, which describe 
the single electron spectrum and the coherence between 
electrons respectively. The occupation numbers of the 
electronic excitations are contained in the Keldysh com- 
ponent G K . 

In general these Green's functions are dependent on 
the time, energy, position and momentum of the parti- 
cle: G — G(E,t,r,p). However typical variations occur 
on a much slower length scale than the Fermi wavelength. 
The Green's functions are sharply peaked around the 
Fermi momentum p = pp, and a considerable simplifica- 
tion can be obtained by integrating G over all momenta. 
A second simplification arises from the short mean free 
path in dirty superconducting films, which averages out 
any dependence on the momentum direction. The result- 
ing equations were obtained by UsadeP-S and they only 
contain what is called the quasiclassical Green's func- 
tions, g{E,x,t) and f(E,x,t). Our experimental obser- 
vations indicate that relevant solutions are stationary, 
so in addition we neglect all time dependences in the 
equations. This choice is partially supported by theo- 
retical work of Snyman et al.^ who demonstrate for a 
simplified system that the solutions for a DC bias are 
always stationary. To parametrize g(E, x) and /(-E, x) 
we use a complex pairing angle 9(E, x) describing corre- 
lation between electrons and holes, and a complex phase 
X- g = cos(6*),/i.2 = sin (6*)e ±ix p21 The normalization 
condition g 2 + /|/ 2 = 1 is automatically fulfilled, while 
the variations of 9{E 1 x) and x(E,x) arc determined by 
the following diffusion equations: 
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with D the normal state diffusion constant. The first 
equation describes how the presence of a local supercon- 
ducting order parameter A(x) generates pair correlations 
9(E,x), which allows to calculate the local density of 
states (DOS) N(E,x) = 5icos(6>). The second equation 
relates the phase gradient of the gap to the presence of 
supercurrents. 

A convenient description of a non-equilibrium super- 
conductor is obtained by introducing a generalized dis- 
tribution function h(E), defined as G K — G R h(E) — 
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h(E)G A . To disentangle the influence of the distribu- 
tion function on the amplitude and the phase of the or- 
der parameter, h(E) is split in the even part (energy 
mode) in particle-hole space /l(E,x), and the odd part 
(charge mode) fr(E,x). The total electron distribu- 
tion functions f(E,x) is then obtained from 2f(E,x) — 
^- — fh{E, x) — /t{E, x). The presence of a charge mode is 
related to the presence of a charge Q* integrated over all 
excitations, and the consequence of inhomogeneity in the 
superconducting system, leading to conversion of quasi- 
particle current to supercurrent. Charge imbalance has 
been studied thoroughly at temperatu res cl ose to Tc, i.e. 
for long lived quasiparticle excitations EEH However, the 
concept of charge imbalance also applies to short-lived 
evanescent states }^ fo r small injection voltages and at 
low temperatures EZHIl 

Conservation of energy E and charge Q result in two 
coupled diffusion equations for and fx- 
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where Hl,t,x are generalized diffusion constants, j e is 
the spectral supercurrent and Rl,t determine the mag- 
nitude of the source term on the right hand side of Eq. 
([3]). The energy current is dominated by the diffusion 
of the energy mode fi . Our Al wires are relatively short 
which means we can neglect inelastic processes™ and the 
energy current is a conserved quantity at every energy. 
The charge current consists partly of a normal current 
driven by a gradient of the charge mode, /„ = IIxV/t, 
and partly of a supercurrent related to a gradient of the 
phase I s = fhjt- Conversion of a normal current into 
a supercurrent implies a change of V/t, and is propor- 
tional to Rt ~ A in Eq. 

The position dependent potential in the superconduc- 
tor e<j)(x) is obtained by integrating the charge of the 
quasiparticle excitations over all energies: 
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In order to conserve charge neutrality, the presence of the 
net charge in the excitations is compensated by a shift 
in the pair chemical potential Sfi cp (x). This means that 
the static electric field E = V^> which drives the normal 
current /„, does not influence the condensate, since it is 
exactly balanced by S/i cp (x) = —e<p(x). If this were not 
the case, the Cooper pairs would accelerate. 



The retarded and kinetic equations ([2] |3| are com- 
pleted with the self consistency relation for A(x): 



A(x) 



NoV f 



Ai 



dE 



sm8e lx - smd*e lx * ) f L - (sm0e lx - smd*e lx * ) / T .(6) 



The charge mode is directly related to the observed po- 
tential drop over the superconductor through Eq. ([5]), 
the energy mode /l only appears implicitly in the gap 
Eq. (§. " 



III. POSSIBLE SOLUTIONS 

In this section we present the numerical solutions of 
Eqs. ([2] - [6j) for the model system shown in Fig. [I] The 
wire can be considered to be one dimensional, as the 
width and thickness are smaller than the dirty supercon- 
ducting coherence length w,t < £ = yj The normal 

equilibrium reservoirs act as boundary conditions, both 
for the superconducting pairing- angles 8 = Vx = and 
the distribution functions /l,t- Temperature enters the 
problem only through the boundary conditions for Jl and 
/t, while all non-equilibrium processes in the wire itself 
are contained in the distribution functions. After an ini- 
tial guess for A(x), the superconducting angles 8 and \ 
are calculated from the retarded equations ph. Subse- 
quently the kinetic equations ([3| can be solved to obtain 
/i and fx- Finally the value of A(x) is updated using 
Eq. ([6]), and this process is repeated until all values con- 
verge. We find two distinct superconducting solutions for 
the problem: (a) one global superconducting state (Fig. 
|2^,) and (b) a bimodal superconducting state separated 
by a normal valley (Fig. [21: 



A. One global superconducting state 

The first solution is characterized by one coherent su- 
perconducting state which extends over the full length of 
the wire, though the strength of the superconducting gap, 
A, is suppressed at the edge of the wire by the presence of 
the normal reservoirs (Fig. [2^,). Though fully supercon- 
ducting, the wire has a finite resistance due to the con- 
version of a normal current into a supercurrent, as shown 
by the position-dependent voltage V. Normal electrons 
which are injected from the metallic reservoirs, decay into 
Cooper pairs over roughly one coherence length £. The 
excess charge Q* associated with the charge mode fx of 
these evanescent quasiparticle states results in the pres- 
ence of an electric field in the superconductor, and hence 
a potential drop over the same length the supercurrent 
increases. These processes correlate with the picture of 
electrons being injected at energies E s» eV, leading to a 
two-step distribution f^, as shown previously by Keizer 
et al. ' While the charge mode non-equilibrium fa relaxes 




(a) One Global Superconducting State (b) Bimodal Superconducting State 
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FIG. 2: (a) The complete wire is in a single superconducting 
state with order parameter A(:r). However near the normal 
reservoirs the condensate carries only a small fraction J s of the 
current as a supercurrent, which results in a resistance and a 
voltage drop at the ends of the wire, over roughly a coherence 
length. At the lowest temperatures a small proximity effect 
can occur at the connection of the bilayer reservoirs to the 
wire (schematically illustrated by dotted black lines), (b) Two 
distinct superconducting domains at the ends of the wire are 
separated by a normal region in the centre of the wire. Due 
to the small supercurrent, the voltage profile is almost equal 
to the normal state. 

over a length scale of £, because of interaction with the 
condensate, the energy mode Jl remains constant over 
the length of the wire due to the absence of inelastic in- 
teractions (Fig. [3^). 

For increasing voltages there is hardly any change in 
the profiles of A, (f>, J s , n , fL,T, until the wire switches to 
the normal state. For example there is no gradual expan- 
sion of the volt age- carrying parts at the end of the wires, 
as one would guess qualitatively. A careful analysis^ in- 
dicates that the energy mode fi, triggers this transition, 
while the current is still far below the critical pair break- 
ing current I c q. 

In performing these numerical calculations we assumed 
that the reservoirs are fully normal down to the lowest 
temperatures. The dashed lines in Fig. [2^ however show 
a schematic picture of a situation where the reservoirs 
are proximitized by the wire, which in fact is a situation 
we encounter in the experiments. The conversion and 
voltage drop occurs primarily in the contact pads, and 
the measured resistance is largely a spreading resistance 
of the contact pad. We will show experimentally that the 
latter contribution can be quenched by the bias and by 
a magnetic field. 

B. Bimodal superconducting state 

A second solution was inspired by our experimental 
results. It consists of two separate superconducting do- 



FIG. 3: The even mode /t and odd mode fr of the non- 
equilibrium distribution function f(E,x). (a) For the global 
superconducting state, a two step distribution is present 
through the full wire, while the charge mode is only present 
at the edges, (b) A strong, non-thermal energy mode non- 
equilibrium /t suppresses superconductivity at the centre of 
the wire. 



mains located at each end of the wire (Fig. |2J)). A strong 
energy mode //, suppresses superconductivity in the mid- 
dle of the wire, while the presence of the cold reservoirs 
near the ends of the wire favors locally the emergence of 
a gap. Modeling this state is complicated, as the pres- 
ence of two superconducting regions gives potentially rise 
to time-dependent processes. We can however avoid this 
complication by assuming that the centre of the wire is 
fully normal. In that case it is possible to proceed nu- 
merically by splitting the wire in two half- wires and treat 
them independently, using 8 = V% = as boundary con- 
ditions. While the distributions at the end of the wire are 
again given by the equilibrium reservoirs, in the middle 
of the wire we match the distribution functions fi and 
fx and their derivatives. The occupation of electronic 
states with energies E + eV, E — eV are coupled by the 
applied voltage, while previously they were independent. 
In addition the superconducting potential mixes parti- 
cle and hole states, and one retrieves relatively complex 
solutions for fL^T(E,x) (Fig. |3b). At the centre of the 
wire the energy mode non-equinbrium is close to a ther- 
mal one, but at an elevated temperature similar to a 
parabolic temperature profile. The remaining structure 
is in essence due to energy-conserving Andreev reflection 
processes .similar to the electron distribution in a SNS 
structured 

The emerging superconducting blobs at the end of the 
wire are relatively small, both in magnitude |A| « |A |/2 
and in size L$ ~ 4£. Due to their limited size only a tiny 
fraction of the total current is converted into a supercur- 
rent, and the voltage profile is almost identical to the nor- 
mal state. While the local microscopic properties at the 
end of the wire show a strong superconducting signature, 
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TABLE I: Overview of the properties of the different samples: L - length, w - width, t - thickness, R n - normal state resistance, 
p - resistivity, D - diffusion constant, £ = ^JhD /2/S. - coherence length, T c - critical temperature, R s - low temperture resistance 
in the superconducting state. For samples indicated with an asterisk, there is no measurement available for T c . We assumed 
the same value for T c as for sample 3 which was fabricated under the same conditions. 



the global properties of the wire are hardly influenced. 
This is true for the current (which is almost completely 
normal) and the voltage profile, but also for the density 
of states and the distribution functions. Apart from some 
small modifications, the distribution function in the wire 
is given by a two-step function. The non-equilibrium en- 
ergy mode Jl is the strongest in the centre of the wire, 
and is the main reason why the superconducting state 
nucleates near the equilibrium reservoirs. The influence 
of /t is limited as the condensate carries almost no (su- 
per)current. 

IV. SAMPLE DESIGN, FABRICATION AND 
CHARACTERIZATION 



purity Al is deposited through a slit in the suspended 
mask to create the superconducting wire and the thin 
bottom layer of the pads. Evaporation of a thick (200 
- 500 nm) copper layers under an angle, which avoids 
deposition through the slit, completes the normal bilay- 
ers forming the reservoirs. Subsequently the Al is oxi- 
dized during 5 minutes in a pure O2 atmosphere with 
a pressure of 4.6 mbar to create an A\O x tunnel bar- 
rier of R n A 300 0,11m 2 . The Cu probes are deposited 
during the last evaporation step under a second angle. 
The size of the wires is measured using scanning electron 
microscopy. The thickness was obtained from a quartz 
crystal monitor used during the deposition of the Al film, 
and calibrated by atomic force microscopy. 



Figure [TJi shows a typical superconducting Al nanowire 
contacted by two massive normal reservoirs, consisting of 
the same thin Al layer covered by a thick Cu layer. For 
reasonably clean interfaces the inverse proximity effect 
of the thick Cu drives the Al normal down to the lowest 
temperatures. The massive volume of the contacts guar- 
antees that they act as equilibrium reservoirs from which 
electrons are injected into the wire. When a bias is ap- 
plied however, the temperature of the reservoirs electron 
distribution function fo(E,T) might for increasing volt- 
age deviate from the bath temperature according toP^ 
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where L is the Lorenz number, Rq the sheet resistance 
of the contact, and r and r\ respectively the electron- 
electron and electron-phonon inelastic mean free path. 
The temperature increase can be considerable, and the 
most obvious way of decreasing it is to minimize the ra- 
tio R\j/R w i re by using thick reservoirs, which we have 
implemented in our sample design. 

The samples are realized by three angle shadow evap- 
oration through a suspended resist mask (PMMA/LOR 
double layer). The parameters of the different samples 
are summarized in Table HI First 50 - 90 nm of 99.999% 



A. Linear response of the nanowire 

Figure [6^, shows a typical current voltage curve (IV) . 
The linear regime extends up to a critical current desig- 
nated by I c \. This initial slope has been measured as a 
function of temperature with an AC technique leading to 
the results shown in Fig. [4] We used a bias current /12 
of 1 jitA modulated at 342 Hz (terminals labels are shown 
in Fig. [T|d). The two point resistance of this 1.4 /im long 
wire (sample la) as a function of temperature displays a 
well-defined pattern (open squares). The spreading re- 
sistance of the contact pads adds a small but finite con- 
tribution of approximately 20 mfi to the measured two 
point resistance. Clearly, at high temperatures the wire 
is normal and has a resistance R n . When the temper- 
ature is decreased below T c — 1.05 K the resistance of 
the wire drops considerably as it becomes superconduct- 
ing. This critical temperature is depressed compared to 
the intrinsic critical temperature of the aluminum due to 
the proximity-effect, as discussed by Boogaard et alPFor 
intermediate temperatures (500-800 mK) the resistance 
appears to saturate at a value Rs ps 1 Q. As we will an- 
alyze further, this is the result of a normal current pen- 
etrating into the wire over roughly one coherence length 
£, yielding a R s ps 2p£/A, with A the crosssection. 

Further lowering of the temperatures leads to a further 
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drop in resistance to almost zero, suggesting that the bi- 
layer contacts are becoming superconducting. To check 
this hypothesis we measured the resistance of identical 
Al/Cu bilayer strips down to the lowest temperatures, 
and find that they stay normal. Instead we attribute 
the vanishing resistance due to the proximity effect by 
the nanowire on the contact point between the normal 
reservoir and the superconducting wire itself. The super- 
conductivity gradually spreads out into the bilayer, lead- 
ing to a normal-superconducting boundary which moves 
from the nanowire into the contact pads (Fig. [2^,). Con- 
sequently, the current conversion resistance in the wire it- 
self becomes gradually less relevant. As the cross-section 
for the conversion moves into the contacts it becomes 
larger, reducing its resistive contribution. Hence, only a 
part of the spreading resistance (~ 20 mf2) is measured. 

This observed pattern changes if we measure the resis- 
tance for a small DC bias current, larger than l c \. Then 
superconductivity in the weakly proximitized region in 
the pads is suppressed, and the differential resistance 
stays constant down to the lowest temperatures (blue 
triangles Fig. |4^i). To further test this hypothesis we 
measured the IV of the wire while we apply a small mag- 
netic field of 7 mT, parallel to the pads but perpendicular 
to the wire (Fig. [Ijj). The vanishing resistance at zero 
bias is no longer observed, while the differential resistance 
at higher biases is identical to the one without magnetic 
field. This indicates that such a small field does not in- 
fluence the properties of the wire, and only quenches the 
weakly proximitized region in the pads. Only at a much 
higher field B ss 100 mT we observe a change in the 
differential resistance of the wire (Fig. |1J;). 

The dashed gray line in Fig. [4] shows the calcu- 
lated two point resistance. The bulk critical tempera- 
ture T cQ = 1.23 K was the only free parameter in the 
fit, while the diffusion constant D = 98 cm 2 /s was ob- 
tained through the relation D = p/Noe 2 . The resistivity 
p is deduced from the normal state resistance R n , using 
7V = 2.2 ■ 10 47 J^m^ 3 for the density of states at the 
Fermi levelP^ The superconducting coherence length is 

obtained from £ = J Although the numerical cal- 
culation agrees quite well with the data (for / > I'd), 
the model overestimates the residual resistance at low 
temperatures. This indicates that the assumption of 
completely normal contact pads is too rigid, as also ob- 
served by Boogaard et alP To include the geometric out- 
diffusion of coherent electrons into the normal pads, we 
adjust the boundary conditions at the ends of the wire 
to: Vf? = —9 /a, which indicates the dilution of super- 
conductivity into the normal pads over a characteristic 
length scale a. With a w 18 nm (full line in Fig. |4| we 
find excellent agreement with the observations. The key 
parameters are listed in Table[j]for the different samples. 
It demonstrates that the linear response of the wires is 
well understood, but that the boundary conditions are a 
sensitive part of the problem even for the thick and wide 
contact pads used. However for bias currents I > I c i the 
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FIG. 4: (a) The two probe resistance versus temperature of 
a 1.4 fim long wire (sample la). Due to the proximity effect 
of the wire on the normal reservoirs, the resistance becomes 
negligible at low temperatures. This weak proximity effect 
can be suppressed by applying a small bias current (b) or 
small magnetic field (c) (sample 4, 200 mK). This 'corrected' 
wire resistance is constant down to the lowest temperatures 
(magenta squares of panel a). A model (dashed line) with 
rigid normal boundary conditions for the pairing angle 9 = 
slightly overestimates the observations. A weaker boundary 
condition (full line), in which 6 decays gradually to zero over 
a characteristic length a shows excellent agreement with the 
experiment. 



system is in a well defined state, which can be connected 
to the theoretical predictions. 



B. Characterization of the tunnel probe 

To measure locally the density of states, the electro- 
static potential ecj>(x) and the chemical potential of the 
condensate p cp we use a normal tunneling probe. The 
current flowing from a normal tunnel probe contacted to 
a non-equilibrium superconductor at a position x is given 
by: 

1 



I T (V,x) = —l m{cos(6(E,x))} 

e -<Vi J— oo 

{f$(E,x)-f$(E + eV)}dE, (9) 



with 6(E, x) the pairing angle, f^(E, x) the charge mode 
non-equilibrium distribution in the superconductor, and 
f£ = 1 - f (E + eV) - f (-E + eV) the distribution 
function in the normal probe. Using Eq. ^ we can 
rewrite this to: 

1 f°° 
I T {x) = —-{ & f>{x)+ I N(E,x)f$dE}. (10) 

€-rt?2 ,/_(Yl 
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The tunnel current consists of two contributions, the first 
one does not depend on the applied voltage, but is com- 
pletely determined by the charge imbalance in the su- 
perconductor, leading to the local electrostatic potential 
4>(x). The second contribution is given by the convo- 
lution of the local DOS N(E,x) of the superconductor 
and the distribution function f^(E + eV) in the normal 
metal. At low temperatures the differential conductance 
of the tunnel contact is a direct measure for the density 
of states in the superconductor. The condensate chemi- 
cal potential of the superconductor can be obtained from 
N(E,x), which is symmetric around E = /i cp . 




0.00 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
Probe Voltage (mV) 

FIG. 5: The differential conductance of the local tunnel probe 
as a function of applied voltage (magenta triangles). Good 
agreement between experiment and theory (full line) is ob- 
tained when the series resistance of the set-up is included 
(blue squares). 



Figure [5] shows a typical measurement of the differen- 
tial resistance for a tunnel probe located at a distance 
of 320 nm = 2.4£ from the normal reservoir of a 4 /im 
long wire. The nanowire is biased just above I c i, to en- 
sure it is in a well defined state. The bias current needed 
to drive the probe is typically four orders of magnitude 
smaller than the bias current of the nanowire, due to the 
high normal resistance of the tunnel junction (Rt — 43 
kf2). Hence it is safe to assume that the properties of 
the nanowire are not influenced by the measurement of 
the probe. One recognizes the coherence peaks at the gap 
voltage, however the subgap DOS is increased in compar- 
ison to the BCS values due to the presence of the normal 
banks, and the driving of the nanowire. The simulated 
local DOS, for the set of parameters, is in good agreement 
with the data, but near the gap voltage a small discrep- 
ancy exists. We attribute this to a series resistance in 
the wiring of the tunnel probe and can correct our data 
for this contribution. We obtain a good agreement be- 
tween the data and the theory using a series resistance 
of Rs = 1.2 kf2, which is the estimated wiring resistance 
of the experimental set-up. 



V. TWO-STATE ANALYSIS AND DISCUSSION 

Wc have realized and studied a total of 7 samples with 
parameters shown in Table |Tj All displayed similar be- 
havior. The non linear current voltage characteristic of a 
typical sample is shown in Fig. [6^,, with two clearly dis- 
tinguished branches. Before discussing the details we first 
indicate the various signatures for processes, which dom- 
inate the various regimes. Increasing the current from 
zero bias we pass I c \ the current at which the proxim- 
itization in the banks is quenched as discussed in Sec- 
tion IV. A. Beyond J cl until 7 c2 , we claim that the wire 
remains in the global superconducting state, character- 
ized by a low and almost constant differential resistance 
R s (Fig. J6Jd) . This resistance reflects the conversion 
of a normal current into a supercurrent and is located 
at the edges of the wire. At the current I = I C 2 the 
wire switches into the normal state, leading to an abrupt 
switch of both the voltage and the differential resistance, 
followed by a constant differential resistance equal to the 
normal state resistance. 

Decreasing the current from the normal state, a kink 
in the measured voltage signals a more subtle transition 
at 7 = 7c3. The measured voltage shows a small deficit 
with respect to its normal state value (black dashed line) 
suggesting the nucleation of superconductivity. We claim 
that superconductivity nucleates here at the ends of the 
wires close to the contact pads in agreement with Fig. (2]d. 
The sudden transition at I = I ci is due to the transition 
from the bimodal to the global superconducting state. 

A first experimental indication to support this inter- 
pretation is provided by Fig. [Hj:, which compares lo- 
cal measurements with measurements over the full wire. 
It shows the two-point resistance of the wire V12/J12 
(squares) as a function of bias current J12, together with 
the apparent resistance V13//12 (triangles) at the ends 
of the wire (see inset for the probe-position and termi- 
nal labels). The probe voltage is multiplied by two for 
comparison, as a similar contribution is present at the 
other edge of the wire. For the lower branch, the as- 
sumed global superconducting state, one observes that 
the voltage drop V13 over the end of the wire is almost 
identical to half of the complete voltage drop over the 
wire-length, a direct proof that this resistance is located 
at the ends of the wire. 

In contrast, in the normal state, the voltages V12 and 
V13 are, as expected, proportional to their respective 
lengths along the wire, V = pL/A. Upon decreasing 
the bias below 7 c3 , where we assume the bimodal state 
exists, one observes over the full length of the wire a 
decreasing resistance for decreasing bias, signaling the 
growing strength of superconductivity somewhere. The 
measured resistance over the end of the wire, however, 
increases compared to the normal state V13//12 > R n - 
Though counter-intuitive this is consistent with the gen- 
eral non-equilibrium present in the superconductor. In 
the following we make a detailed analysis of both super- 
conducting states, and place the experimental results in 




FIG. 6: The voltage V12 (a) and differential resistance (b) of a 4 /im long wire (sample 3b) as a function of bias current J12, 
measured at 200 mK. We define four different regimes with boundaries labeled I c i — Id, each characterized by a nearly constant 
differential resistance. The critical currents I C 2 and Id are defined as the bias currents where the wire switches between the 
two hysteretic voltage branches. I c i and I c i are the transition points between the two different states of one branch, (c) 
The apparent resistance of the complete wire V12//12, and of the edge of the wire V13//12 as measured with a voltage probe, 
multiplied by two for the ease of comparison. 



the context of the theoretical model. 



A. Global superconducting state 

Figure [7] shows two-point measurements of the lower 
branch of a 1.4 ^m long nanowire (Sample la) at three 
different bath temperatures. In view of the analysis 
shown Section I IV Al we assume that the resistance of 
the wire is primarily determined by the charge mode 
of the distribution function fx(E,x), which depends on 
the position dependent density of states and the order 
parameter A(x). The weak dependence of the differen- 
tial resistance on the current indicates that the super- 
conducting properties of the wire hardly change with in- 
creasing bias (open symbols). Though numerical sim- 
ulations (filled symbols Fig. [7| show the same qual- 
itative behavior, the simulations seem to overestimate 
the bias current at which the differential resistance be- 
gins to increase. Hence the observed switching current 
I C 2 is also slightly lower than predicted. Ignoring this 
small discrepancy, the simulated data show good agree- 
ment with the experiment over the complete temperature 
range (inset Fig. [7|. At the same time the observed val- 
ues for the critical current (or critical voltage) are much 
smaller than what one would expect for a pair-bre aking 
current, experiment allypHMl as we u as theoretically! 25 * 26 ^ 
This demonstrates that the non-equilibrium processes 
should be taken into account in evaluating the param- 
eters. The remaining deviations between theory and ex- 
periment suggest, most likely, that the temperature of the 
reservoirs deviate from the bath temperature for higher 
driving currents, as expected from Eq. ([Sj). 

Figure [8^1 shows the differential conductance dl^/dV^ 
of a tunnel probe, located at a distance of 320 nm (=2.4 
£) from the normal reservoir of a 4 ^m long wire (Sam- 
ple 3b). The wire is biased at a fixed current I12 with a 
corresponding voltage Vx%. At the same time the probe 
current /13 is varied while measuring the probe voltage 
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FIG. 7: Two probe voltage (a) and differential resistance 
(b) as a function of bias current, for a 1.4 (im wire, at three 
different bath temperatures. Open symbols: experimental 
data. Filled symbols: numerical simulations. The critical 
current as a function of temperature (inset). 



V13 (Fig. [TJd). The evolution of the local DOS for in- 
creasing bias is shown on the right side in Fig. [8^,. The 
conductance at zero bias V p = increases slightly for 
increasing bias, while the coherence peaks get further 
smeared out. The dependences are however weak, and 
even right before the switching current I c2 the DOS is 
hardly affected by the drive current. These observations 
are in close agreement with the theoretical predictions 
and confirm the idea that the superconducting state re- 
mains globally stable. For increasing bias the resistance 
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remains located at the ends of the wire and the DOS does 
not change either. Unfortunately, we have not been able 
to directly measure Nevertheless we believe that 

this energy-mode non-equilibrium triggers the transition 
at I C 2 as analyzed by Keizer et alP 



B. Bimodal superconducting state 
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FIG. 8: The local density of states for (a) the global supercon- 
ducting state and (b) the bimodal state, for different bias cur- 
rents 7i2 of the nanowire, measured at 200 mK. For the global 
superconducting state the gap is only weakly dependent on 
the bias current, while for the bimodal state one observes a 
DOS gradually changing from a normal into a superconduct- 
ing state. 

The continuous transition, with decreasing bias, from 
the normal into a superconducting state at (Fig. [6]) 
indicates that the emerging superconducting state is ini- 
tially very close to the normal state. For lower bias cur- 
rents the absolute resistance gradually decreases (Fig. 
[6j;) , which suggests that an increasing fraction of the cur- 
rent is carried by the emerging condensate (Fig. [2Jd). A 
similar picture is observed for the local density of states, 
plotted in Fig. [8}d for different bias currents I\2 of the 
nanowire. Below I C 3 a gradually increasing gap is found, 



unambiguously showing the emergence of superconduct- 
ing order. Close to J C 3 the DOS at the position of the 
probe evolves in a continuous way from a flat spectrum 
into a spectrum with a gap. However at J C 4 one ob- 
serves an abrupt transition to a situation with a stronger 
gap. The abruptness indicates that it is a transition 
from two distinct superconducting states, which directly 
proofs that at least two microscopically distinct super- 
conducting states exist. Though the simulations for the 
local DOS agree well at currents close to 7 C 4, they do 
not account in detail for the gradual evolution between 
the normal and superconducting state at I c ^. At this 
point we assume that the reservoirs start to heat up, 
and can no longer be treated as equilibrium reservoirs 
with T — Tfr. Overall, the model supports the picture of 
the emergence of the superconducting state quite nicely, 
with the strongest non-equilibrium in the wire occurring 
at 7 C 4, with the reservoirs most closely to equilibrium at 
T = T b . 

The electro-chemical potential of the superconducting 
condensate, fi cp , is determined from the minimum of the 
measured DOS /i cp = eV\ m i n (£>os)- It is found that at 
this probe-position fi cp is equal to the electrostatic po- 
tential Vi of the adjacent reservoir. Measurements with 
a probe in the middle of the wire show that, in the same 
bias regime, the voltage is equal for both sides of the 
wire, which means that the state is symmetric, and a sim- 
ilar superconducting region should exists near the other 
reservoir at a potential Vi. If these two regions were 
part of one global superconducting state, there would be 
a voltage drop AV — V\ — Vi over the superconduct- 
ing potential /i cp of this state, and a superconducting 
phase-slip process should occur. However, according to 
the Josephson relation 2eV = d\/dt it would be at a 
frequency v w 8H/A., which is too high compared to the 
energy-gap. 

The fact that the two-point resistance is so close to 
R n , the gradual increase of the DOS at the position of 
the probe, and the electro-chemical potential of the con- 
densate demonstrate that two separate superconducting 
regions emerge at the edges of the wire. The physical rea- 
son is the energy mode non-equilibrium, as discussed by 
Keizer et aP for the lower branch, but similarly for this 
upper branch. At the bias I C 4 the wire is still largely nor- 
mal and /e is given by the two step distribution-function. 
In the middle of the wire the width of the step is several 
times bigger than the superconducting gap. Through re- 
lation Eq. ^ it is seen that this suppresses fully the 
nucleation of a gap, while the cold equilibrium reservoirs 
favor the emergence of a gap at the edges of the wire. 
Simply put, the ends are cold where the centre of the 
wire is hot. Therefore we conclude that the results are 
most easily understood as due to two distinct supercon- 
ducting domains, separated by a normal central region, 
what we have called the bimodal state. 

Finally, we discuss the voltage V13 measured by the 
probe when the wire is biased into the bimodal state (tri- 
angles Fig. |6|;) . Close to T c the voltage measured by such 
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a normal probe is equal to the electrochemical potential 
of the quasiparticle bathP^ At low temperature^^] anc l 
for short-lived quasiparticles^ it is impossible to define 
a quasiparticle bath with a well-defined chemical poten- 
tial, however the measured voltage is still related to the 
local electrostatic potential ecf) (using Eq. 10): 



N(E)0(E + eV)dE. 



(11) 



For a relatively small charge imbalance ecj)(x) <C A, the 
measured voltage equals the local electrostatic potential 
e(f>(x) divided by the local DOS in the superconductor 
at zero energy: V ~ cf>/N(0). Hence the voltage mea- 
sured with the tunnel probe can be larger than the local 
potential e<p(x). 



VI. CONCLUSION 

We have analyzed a well-defined model-system of a 
superconducting wire between two massive normal con- 
tact pads. We demonstrate that this NSN system, when 
driven by a current, has two distinct mctastable super- 
conducting states. 

For low bias we find a global superconducting state 
with most of the resistance occurring as a current- 
conversion resistance at the ends of the superconducting 
wires where normal current enters. Although resistive we 
demonstrate that the whole wire including the edges con- 



tinues to be in one coherent superconducting state. This 
state does hardly change for increasing current, until the 
wire switches abruptly to the full normal state at a cur- 
rent, which is much lower than the critical pair-breaking 
current. On a microscopic level the distribution function 
changes considerably and is strongly different from the 
commonly used parabolic temperature profile. A numer- 
ical analysis based on the non-equilibrium quasiclassical 
Green's functions shows that the switching current is de- 
termined by the non-equilibrium electron distributions, 
in good agreement with the experimental results. 

For high bias, decreasing the current from a fully nor- 
mal state, we find that the superconducting state emerges 
as two decoupled domains at the ends of the wire. The 
vicinity of the cool equilibrium reservoirs favors the nu- 
cleation of the superconducting state at these ends, while 
strong non-equilibrium at the center of the wire continues 
to suppress the superconductivity. Upon further lower- 
ing of the bias current the two domains grow in strength 
until the wire switches back to the low resistive, globally 
superconducting state. We speculate, that the transition 
from one state to the other, is triggered by a condition 
in which the Josephson coupling energy between the two 
domains exceeds the thermal energy at that bias point. 

We would like to acknowledge Nanofridge (Nanosci- 
ERA), Microkelvin (No. 228464, Capacities Specific Pro- 
gramme) and the Foundation for fundamental research 
on matter (FOM). 
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